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We use 2-color QCD as a model to study the effects of simultaneous presence 
- - of the so-called 9 parameter, chemical potentials for baryon number, /x^ and for 

Q ' isospin charge, /i/. We pay special attention to 6, fj^B, fJ'i dependence of different 

O ■ vacuum condensates, including chiral and diquark condensates, as well as the gluon 

^ condensate, (^^G^^C^*^"), and the topological susceptibility. We find that two 

phase transitions of the second order will occur when 9 relaxes from = 27r to = 0, 
if n is of order of the pion mass, m^. We demonstrate that the transition to the 
superfluid phase at 9 = it occurs at a much lower chemical potential than at 9 = 0. 
We also show that the strong 9 dependence present near = vr in vacuum (Dashen's 
phenomenon), becomes smoothed out in the superfluid phase. Finally, we comment 
'sj- ] on the relevance of this study for the real world with Nc = 3. 

o 

^ , PACS numbers: 

O 

o 

^ ; I. MOTIVATION 

Ph. 

In this paper we investigate the behavior of 2-color QCD under the influence of three 
^ . parameters: 9, fis and /x/. The main motivation for such a study is, of course, the attempt 
to understand the cosmological phase transition when 6, being non-zero and large at the 
very beginning of the phase transition, slowly relaxes to zero, as the axion resolution of 
the strong CP problem suggests. Therefore, the universe may undergo many QCD phase 
transitions when 6 relaxes to zero. Another motivation is the attempt to understand the 
complicated phase diagram of QCD as a function of external parameters 6, fiB and fij. 
Finally, our study may be of interests for the lattice community - the determinant of the 
Dirac operator for = 2 is real when ^ = tt in the presence of nonzero /i. As we show, in 
this case the superfluid phase is realized at a much lower chemical potential than at ^ = 0. 
This gives a unique chance to study the superfluid phase on the lattice at a much smaller /x 
than would normally be required. 

To study all these problems in real 3 color QCD at finite fis is, of course, a very difficult 
task. To get some insight into what might happen we shall use a controlled analytical 
method to study these questions in the non-physical (but nevertheless, very suggestive) 
Nc = 2 theory. We use the chiral effective Lagrangian approach to attack the problem. We 
shall determine the phase diagram in the fiB,fJ^i,(^ planes, various condensates and lowest 
lying excitations. We expect that our approach is valid as long as all external parameters 
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^B,fJ'i and the quark mass, rrig, are much smaller than Aqcd- We perform most of our 
calculations for the case of two flavors Nf = 2 where the algebra simplifies considerably. 

One exciting effect that we find is that 6 dependence of the theory at fixed fi may become 
non-analytic. This is due to the fact that the critical chemical potential for transition to 
the superfluid phase varies with 9. Therefore, a change of 6 might trigger a second-order 
phase transition, accompanied by a discontinuity in the topological susceptibility x- We 
also find that the strong 6 dependence, present near 6' = vr in vacuum, is washed out in 
the superfluid phase. We expect that for equal quark masses a first order phase transition 
(Dashen's phenomenon) will occur in the N^. = 2, Nj = 2 theory at 9 = n, in the normal 
phase, but will disappear in the superfluid phase. 

We also find some interesting results, which appear even at 6' = 0. Most importantly we 
compute the dependence of the gluon condensate, {^2^G'^^^G^'^'^) , on the chemical potential. 
The gluon condensate decreases with density near the normal to superfluid phase transition, 
but, counter-intuitively, increases for <^ <^ Aqcd- 

We also evaluate novel vacuum expectation values which appear in the superfluid phase: 
(m"^7oC75r2(i) in the baryon breaking phase and {iu'jo'y^d) in the isospin breaking phase. 
These densities, being nonzero even at 6' = 0, nonetheless have never been discussed in the 
literature previously. These densities, themselves, break the baryon and isospin symmetries 
respectively, and so may be considered as additional order parameters. 

The presentation of our results is organized as follows. In section II, we introduce our 
notations for the low energy effective Lagrangian. In section III, we introduce the 6 parame- 
ter into the effective Lagrangian description. In section IV, we discuss the phase diagram of 
our theory in detail, computing the spectrum of lowest lying excitations, characterizing the 
phases in terms of chiral condensates and densities and paying special attention to physics 
near the point 6 = ir. In section V, we check that our results satisfy known Ward Identities 
supporting the self consistency of our approach. In section VI, we study the gluon conden- 
sate {-^^G'Ij^^G'^'"^) as a function of /i and 6. In Conclusion, we discuss the relevance of the 
obtained results for 3 color QCD and make some speculative remarks on evolution of the 
early universe during the QCD phase transition. In appendix I, we clarify some technical 
issues associated with global aspects of the goldstone manifold. 



II. THE EFFECTIVE THEORY AT FINITE fiB AND m 



Two color QCD at zero chemical potential is invariant under S\J(2Nf) rotations in the 
chiral limit. This enhanced symmetry (as compared to the SU(A^/)xSU(A^/)xU(l) of three 
color QCD) is manifest in the Lagrangian if we choose to represent it in a basis of quarks ip 
and conjugate quarks [J, 0]. For Nf = 2 we use. 
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where the Pauli matrices T2 and (J2 act in colour and spin space respectively. We work in 
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Euclidean space and use the definitions, 7,, 







75 



1 



3 



The microscopic Lagrangian then reads, 

L = i^^a^D^^ (2) 

and possesses a symmetry, 

^ ^ f/^, U e SU{A) (3) 

The enhanced symmetry manifests itself in the low energy effective theory through the 
manifold of goldstone modes associated with spontaneous breaking of chiral symmetry, 
SU{2Nf) Sp{2Nf). In our case, Nf = 2, and the goldstone manifold is SU(4)/Sp(4), 
corresponding to the condensation of — SU(4) flavor sextet. The fields on this man- 
ifold can be represented by a 4x4 antisymmetric unitary matrix E, with det S = 1, that 
transforms under Q as, 

S ^ f/Sf/^ (4) 
We parameterize the vacuum manifold as, 

E = f/Sef/^, U e SU{4), K=(^1'q^ (5) 

In what follows we use notations suggested in 0, 0, 1^ for the description of baryonic as 
well as isospin chemical potentials. In these notations the baryon charge of the quark is 1/2, 
which comes from 1/Nc, so that the baryon (diquark in Nc = 2) has baryon charge 1. Thus, 
chemical potentials enter the microscopic Lagrangian as, 

L = ip'y^D^ip - ^fiBi'l^'ip - ^fiiip'y^a^ip (6) 

In the basis of SU(4) spinors ffl the baryon and isospin (third component) charge matrices 
in block-diagonal form are 0,0,1 



so that the Lagrangian reads, 

L = i¥a,D,^ - ^^fiBB + fiil)^ (8) 

The effective Lagrangian for the field S of goldstone modes is determined by the sym- 
metries inherited from the microscopic two-color QCD Lagrangian. To lowest order in 
derivatives and at zero quark mass the effective Lagrangian isjH, 

£ = — TrV.SV.St (9) 

The //-dependence enters the effective Lagrangian through the covariant extension of the 
derivative, 

9oSt ^ VoSt = 9oSt + [{^iBB + + S(/iB5 + = (10) 

required by an extended local gauge symmetry [l[. Therefore, to this order in chiral per- 
turbation theory, the Lagrangian at finite fi does not require any extra phenomenological 
parameters beyond the pion decay constant, F. This fact gives predictive power to chiral 
perturbation theory at finite yU. In using the effective Lagrangian constructed above we 
must, of course, assume that chiral symmetry for Nc = Nf = 2 QCD is spontaneously bro- 
ken. Since we have regarded the hadronic modes as heavy, the theory is expected to be valid 
only up to the mass of the lightest non-goldstone hadron. 
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III. THE MASS TERM AND 6 PARAMETER 

The mass term in the fundamental Lagrangian is defined as, 

Lm = niuUU + m^dd (11) 
while the 9 term in the fundamental Lagrangian is, 

Le = ^9■ ^ (12) 

We keep m„ ^ rud on purpose: as is known = is a very singular limit when one 
discusses 9 dependence, see below. We would like to incorporate the 9 dependence directly 
into the mass matrix. This can be achieved by performing a chiral rotation. 

With this field redefinition, the topological 9 term in the Lagrangian disappears, due to the 
axial anomaly, and the mass term becomes, 

Lm = ij^^^^M^ + ^l^Mi; (14) 



M = e-*^/^/ ( "r^" ° ) (15) 



where the mass matrix M is, 



rrid 

In the basis of SU(4) spinors ([Q), the mass term becomes, 

Lm = ^^^A^aaTs* + h.c. (16) 

where, in block-diagonal form, 

^=[-M ) (1^) 

The transformation properties of Lm under (jS)) imply that to lowest order, Ai enters the 
effective Lagrangian as, 

Cm = -gReTi{MJ:), (18) 
where the coefficient g is determined by the chiral condensate in the limit m ^ 0^ , 9 = 0, 



.-'^ in) 



as will be confirmed below. In our notations the chiral condensate includes the sum over all 
fiavors, (?/'^) = Y.f {i'fi^f)- 

The chiral effective Lagrangian incorporating the effects of fiB,fJ'i,9 and non-zero quark 
masses, thus becomes, 

£ = — TrV^SV^S^ - ^i?eTr (ME) (20) 
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We shall use the Lagrangian ()2Up for the rest of this work. 

So far our discussion easily generalizes to arbitrary Nf. However, significant algebraic 
simplification can be obtained by considering Nf = 2. Indeed, for Nf = 2, the effective 
Lagrangian with m„ 7^ m^, 0, can be reduced to the same Lagrangian but with 
= m'^ and 6' = 0. This is achieved, by performing an SU{A) (more specifically SU{2)a) 
rotation, 

S = f/oSf/o^ (21) 

with the particular choice of, 

{rriu + rrid) cos(6'/2) 

cos a — 



sma 



v/(m„ + md)2cos2(^/2) + {niu - mrf)2sin {6/2) 
{rriu - rrid) sm{6/2) 



v/(m„ + md)2cos2(0/2) + 
Our parameter a is related to the commonly used Witten's variables 0„, 0^0, via. 



(22) 



0„ = ^/2 - «, 0, = 9/2 + « (23) 

+ 0d = 9, TTlu sin 0„ = TTld siu 0rf (24) 



After such a transformation, the Lagrangian (pn|) takes the form. 



C = — TrV^SV^S"^ - ^m(^)i?eTr (^Motl^ (25) 



A<o = ( _° n)' "^(^) = ^((m. + mrf)2cos2(^^/2) + (m„-m,)2sin2(e/2))^ (26) 



with. 



1 y ' 2 

The detailed explanation of this reduction, which along the way clarifies certain global prop- 
erties of the vacuum manifold, is presented in apppendix I. Technically, the simplification is 
due to pseudo-reality of SU{Nf = 2) (see also section IVC for a more quantitative discus- 
sion). 



IV. PHASE DIAGRAM 

A. Vacuum Alignment and Spectrum 

Our next step is to find the classical minimum of the effective Lagrangian (j^UI) to deter- 
mine the phase diagram, pattern of spontaneous symmetry breaking and, subsequently, the 
spectrum of excitations. For arbitrary Nf, quark masses, 9 and chemical potentials this is a 
non-trivial algebraic problem. However, as was shown in the previous section, for Nf = 2, 
the effective Lagrangian reduces to the form (j25|) . which was already analyzed in Thus, 
we may immediately read off all quantities of interest. 
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First, let's study the phase diagram for fixed 712^,711^, 9. To get acquainted with our 
theory, let's begin with the trivial environment /i^ = /i/ = 0. The effective Lagrangian 
possesses an S'p(4) symmetry at this point. The classical minimum is given by, 

(S) = (27) 

The 5*^(4) symmetry is unbroken. The low-lying excitations are a quintet of pseudo- 
goldstones (3 pions and 2 diquarks), with dispersions, 

E = ^v^ + ml{e) (28) 

T^AO) - (29) 

The pseudo-Goldstone mass m.,^ acquires a dependence on 9 through the effective quark mass 
parameter m{9) (j26|) (this 9 dependence is implicitly implied in all formulas below, unless 
otherwise stated). As we shall see, the whole phase diagram turns out to be determined by 
the parameter ■m.^{9). We note that m.^{9) reaches its maximum at ^ = and minimum at 
9 = IT. Moreover, for m„ = m^, 9 = tt, m.^ vanishes to first order in M. 

We note that strong P and CP symmetries are explicitly broken in the system with 
9^0. So at 6* 7^ 0, the pions (diquarks) are no longer pure pseudoscalars (scalars). This 
will become particularly clear when we discuss Bose-condensates of our goldstones in the 
superfiuid phase. 

Now let's turn on chemical potentials. For /ib 0, f^i ^ 0, the symmetry of the problem 
is broken to U{1)b x ^7(1)/.^ We introduce the following notations to describe vacuum 
alignment of S at finite chemical potentials. 

As is known[i'], there are 3 distinct phases in the (/x^,/^/) plane, ^ 
I. Normal Phase (N): IfJ^s] < fnn{9), < m.„{9) 

(S) = (31) 
Symmetry breaking: U{1)b x U{l)i U{1)b x U{l)i 

Spectrum: 

E = + ml±fiB 

tt" E = + ml 



TT* E = a/p^ + ± fjj 



(32) 



^ If only one of the chemical potentials is turned on, say /i/ — 0, fiB 7^ 0, then the symmetry is actually, 
SU{2)v X U{1)b. 

^ In the original paper [Sj, a certain physically reasonable ansatz was taken for the classical static minimum 
(S) of ((201 ■ It was shown that this ansatz is, indeed, a local minimum, and the authors assumed that 
this minimum is also global. We note, that using the explicit parametrization of the vacuum manifold 
presented in Appendix I, it is possible to prove that the ansatz is, indeed, a global minimum. 
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II. Baryon Phase (B): l/x^l > mT,{6), < \hb\ 



1 / 4 

ml „ / m 



2 



(S) = -fS,+ 1--/) (33) 



Symmetry breaking: V{\)b x f/(l)/ — > ^(1)/ 

Spectrum: 



t = P^ + ^/i|(l + 3^)±/.j4p^^ + i/.|fl + 3^'''' 

2 V t^B/ \ f^B ^ \ t^B 



(34) 



III. Isospin Phase (I): > ^2^,(9), \^b\ < l/W/l 

(S) = ^S.+ (i-^)'e, (35) 

Symmetry breaking: U{1)b x ?7(1)/ U{1)b 

Spectrum: Same as for B Phase, but with q"^ tt"^ and fis ^ fii- 

The phase transition between N phase and B phase, as well as phase and / phase 
is second order, whereas the phase transition between B phase and / phase is first order. 
As noted in |5i], the symmetry of the phase diagram/spectrum, with respect to yU/ 
is a direct consequence of the symmetry of the microscopic theory, ^ —d^f, dfi ^ di, 

12b <^ /i/. 

Thus, the phase diagram in the (/i^, fij) plane looks the same at 7^ as at 6' = 0, with 
the important replacement, — > m1{6). This is a very natural conclusion. Indeed, at ^ 7^ 
diquarks (pions) still carry baryon (isospin) number. Hence, their energy is lowered at finite 
baryon (isospin) chemical potential. As soon as /i^ (/i/) reaches the vacuum diquark (pion) 
mass m^(6'), Bose-condensation occurs leading to spontaneous breaking of U{1)b {U{1)i) 
symmetry. 

Quantitatively, the 6 dependence of the Goldstone mass mT^{6) implies that the transition 
to superfluid phase is shifted to a smaller chemical potential /ib, /i/, compared to = 0. 
In the limiting case, when iriu = md and 6 = it, the transition occurs in the vicinity of 

/i = (see Section IVC for a more precise discussion). For physical values, = 7MeV, 

1 

ruu = 4:MeV, the transition at 6* = tt occurs at /i = (^ ^^"^^ j ' m^{0) ~ 70MeV. 

B. Chiral Condensates and Densities 

In section IVA, we have established the phase diagram of = 2, Nf = 2 QCD at 
finite Hb-, iJ^i and 9. In this section, we wish to characterize this phase diagram in terms of 
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chiral condensates and densities. Similar computations have been performed^, 1^ at ^ = 0. 
However, we evaluate a wider range of expectation values and find some condensates that are 
non-zero even at ^ = 0, which have not been discussed in the original papers. As expected, 
we also find new condensates at ^ 7^ 0. 

We follow the standard procedure for computing microscopic condensates from the effec- 
tive Lagrangian. We start from a slightly generalized version of the microscopic Lagrangian 
(jHl) together with the mass term (fT^ . 

L = i^V^D^^ - ^^T^ + ^ (^^JasTs^ - ^WaTa^*) (36) 

Here the hermitian, traceless matrix T incorporates the chemical potentials for all 15 
charges associated with the S'f/(4) symmetry, and the chiral condensate source J is an 
arbitrary, antisymmetric matrix (12 real components). We may express T and J in terms 
of a basis, 

T = tAXA, J = JaXa, Xa = ^Xa) = 0, Xj = -X,, tA,Ja G ^ (37) 

Differentiating the vacuum free energy density JF, we obtain our condensates and charge 
densities: 



— = -(^^X,a2r2^-^txt„a2r2^*) (38) 

The relations ()38p hold for any T, J, however, we will apply them when the derivatives and 
expectation values are evaluated at physical parameters, T = Tq = /ib-B + /i// and J = Ai. 

In the effective theory, the sources T and J are incorporated by replacing Tq ^ T in the 
covariant derivative (fTUI). and J\4 ^ J in the mass term (fTS|) . The condensates (|3H|). thus, 
become, 

— = F2(Tr(StA^VoS- VoStA^S)) 

OlA 

— = -(7(i?eTr(X,S)) (39) 

It remains to evaluate the expressions (j39|l . with S given by the time- independent, clas- 
sical minimum of the effective Lagrangian (j2DI). We must remember that to simplify algebra 
we expressed, S = UoT,Uo^ , with Uq given by (j22I)- We should also remember that we in- 
corporated 6 dependence into the mass-matrix by a chiral rotation (fT!?|) of the quark fields. 
The condensates and densities, expressed in terms of the original quark fields, are listed in 
Tables 1,2. We define the charge conjugation matrix C = 707275- We also introduce the 
parameter X{9) in Tables 1,2, 

1 Normal Phase 
X(9) = <j Baryon Phase (40) 



Isospin Phase 



which obtains its 9 dependence through m^{9). 
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TABLE I: Chiral condensates in Nc = Nf = 2 QCD at finite 9 



Condensate/ (V''0)o 


N Phase (A = 1) 


B Phase (A = 

Mb 


I Phase (A - 







-lcos(|)(l-A2)^ 





U^CT2d 





-isin(f)(l-A2)^ 













-icos(|)(l-A2)^ 


ud 








-isin(f)(l-A2)^ 


uu 


^Acos(| — a) 




— |Asin(| — a) 


dd 


|Acos(| + a) 


id'y^d 


-^Asin(|+a) 



TABLE IL Densities in Nc = Nf = 2 QCD at finite 6 



Density 


N Phase (A = 1) 


B Phase (A = 


I Phase (A = ^) 







4FVb(1 - A2) 





iu^^oCj5T2d 





-4F2^bA(1- A2)^cos(a) 










-4FVBA(1-A2)5sin(a) 













4FV/ (1 - A2) 


iM7o75(i 








4FV7A(1- A2)^cos(a) 


u'jod 








4FV/A(l-A2)^sin(a) 













(FjoCT2d 
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We can now see, how our phase diagram is described in terms of condensates and charge 
densities. First, let's gauge our intuition by considering the Normal phase. At ^ = 0, the 
parameter a of eq. ()22|) is 0, and the only condensates (Table 1) are (mm) = (rfrf). At non-zero 
6*, we also get condensates {iu'-y^u), (id'y^d), while (uu), (dd) get depleted. The appearance 
of P and CP odd condensates {iu-y^u), {id'~^^d) is a direct consequence of explicit P and 
CP breaking by the 9 term. Finally, for 6' 7^ 0, 7^ m^, the parameter a 7^ 0, and we 
see explicit effects of isospin symmetry breaking: {uu) 7^ {dd) (correspondingly for P odd 
condensates). Such effects are absent to lowest order in M at ^ = 0. As expected, all charge 
densities (Table 2) in the Normal phase vanish. 

Let's now see what happens in superfluid phases. At ^ = 0, the Baryon phase is charac- 
terized by a scalar diquark condensate {iu'^C'j5T2d) , which breaks the U{1)b symmetry. The 
Isospin phase is characterized by a pseudo-scalar pion condensate {iu'j^d), which breaks the 
U{l)i symmetry. These condensates appear at the expense of depleting {iIjiI>). As expected, 
at finite 9, U{1)b and U{l)j violating condensates of opposite parity also appear: {u^CT2d) 
in B phase and {ud) in / phase. 

The Baryon and Isospin phases also carry non-vanishing SU{A) charge densities. The I 
phase, is characterized by the isospin density, 

ni = ^(^70^3^) = (1 - ^) (41) 

This is precisely the density, which one expects to induce by applying an isospin chemical 
potential fij. At ^ = it coincides with the previous results addition, we also 

obtain the following axial charge density, 

UA = {lupoid) = AF'fii'^ (l - ' cosa(^) (42) 

which has not been discussed previously in the literature even at ^ = 0. This is the axial 
charge density, corresponding to off-diagonal generators of the SU{2)a group, which is both 
spontaneously and explicitly broken. Note that the axial charge density ()42j) does not vanish 
already at 6' = 0. Thus, we for now concentrate on ^ = 0, and hence a{9 = 0) = 0, to better 
understand the physical nature of this new density ()42j). For simplicity, we take \^b\ < 

The density ha spontaneously breaks the U{l)i symmetry and, hence, may be considered 
as an order parameter alongside the pion condensate, 

1 - / m^\^ 
{n-) = {^u^,d) = --{^P^P),^^l--^j (43) 

There was no explicit chemical potential conjugate to ua in the Lagrangian - once U{l)i 
is already spontaneously broken by (vr), ha is induced automatically. The quantitative 
behaviour of these two order parameters is somewhat different. The pion condensate mono- 
tonically increases with /// after the Normal to Isospin phase transition, and (vr) — —^{ipip)^ 
for Hi S> mTT. On the other hand, the new charge density ua first increases after the phase 
transition, reaches a peak at yU/ = S^/'^m^, and then decreases to for /i ^ itLt,. Of course, 
we always consider only fii,fiB -C Aqcd- 

One can understand the appearance of a new condensate ha = {iu'-^ol^d) in the following 
simple way. We are in the phase where the isospin density, ni ~ {u'jou) — {d'jod), as well 
as the condensate, {n~) ~ {iify^d), do not vanish. This implies that our ground state can 
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be understood as a coherent superposition of an infinitely large number of vr" mesons. We 
expect that we do not disturb the ground state of the system by adding one of these n~ 
mesons. On the other hand, we can relate the matrix element with an extra vr" meson 
to the matrix element without the vr" using the standard PCAC technique, (A|0|i?7r) ~ 
z(74|[0, Q^]\B). The coefficient of proportionality would not be precisely 1/F in the present 
case because our pions are not in a trivial vacuum, but rather in the (vr") condensed phase. 
However, we expect that the general algebraic structure of the vacuum expectation value will 
be obtained correctly using this approach. Indeed, taking O = U'^qu — d'jod, as the isospin 
density and calculating the commutator [0,(5^], where = J d^xU'jo'y^d = J d^xu^'~^^d is 
the axial charge, one obtains the structure (iu'-fo'-f^d) entering the eq. (021) • Therefore, if we 
expect the vacuum expectation value of (O) = {u'you — djod) to be nonzero, we should also 
expect a nonzero value for the axial density {iu'yo'y^d). This logic is definitely supported by 
the explicit calculations (j42|) . 

One can also test formula (j42p at small isospin density rij. In this case our system may 
be understood as a dilute Bose-Condensate of non-relativistic tt~ particles How is ha 
manifested in this terminology? We shall work at fixed isospin number density (instead of 
at fixed III). Moreover, we will temporarily work in Minkowski space. In the Isospin phase, 
the diquarks are not important as we saw, so we parameterize S as. 



-U 




U E SU{2) (44) 



The field U transforms as U ^ LUB) under SU{2)l x SU{2)r and the effective Lagrangian 
for U reads, 

C = F^ [Tid^Ud^'U^ + 2mlReTiU) (45) 

Thus, we see that the Lagrangian describing the pion sector of N^. = Nf = 2 QCD is exactly 
the same as the one describing = 3, Nf = 2 QCD. We express, U = exp (^^^^)- Similarly 
to eq. (jnHI),®, we identify, 

^PJ^'^i; = 2iF^Ti (^[5^f/, f/^] ~ -e'^^'^^'^vrV (46) 
^T^yV" = -2iF2Tr |^{a^f/,f/t}y^ ^ 2Fd''n'' (47) 
#7'y^ = ^^?Tr (^{U - f/t)y^ ^ -g^ (48) 

where we have expanded the corresponding currents to leading order in tt fields. We also 
expand the Lagrangian to fourth order in vr fields, 

C = -d.nd^'n - -mlTX^ L^vr^S.TrS^vf + -^dJ^^)d^{^^) + -^mliii^f (49) 

2 A* 2 ^ 2AF^ / \ / ' gg^2 ttv / \ i 

We can ignore the 7r° particles as they are irrelevant for vr^ condensation. It is useful to 
combine vr = ^^(vr^ + iTr^). To describe non-relativistic physics involving vr" particles, we can 
replace OqTI imT^ir, diir 0, in the quartic terms of Lagrangian (|49p . Finally, we adopt 
a non-relativistic normalization of our vr" field, by introducing, a canonical, non-relativistic 
Bose field (of dimension 3/2), 

(p^ = V^TT (50) 
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The Hamiltonian density in terms of the field reads, 



H = i^M'^M + m.0t0 + (51) 



while the condensates and densities become 

1 



ni = -tpast/j = (jy'ct) (52) 



Ha = ^7075^ = —2F^/rn^(jr (53) 
iu'^^d = - ^ 0"^ (54) 

In this language, we see that both of our f/(l)/ order parameters, ua and iu'y^d are expressed 
in terms of the same non-relativistic Bose field 0^. 

The energy density of a spatially uniform Bose- Condensate as a function of isospin density 

is, 

e = m^ni + ^^nj (55) 

Therefore, the isospin chemical potential, 

de 1 
= ^ = + Y^-/ (56) 

One can check that agrees to first order in rij with the result (jlT|l obtained in the 
grand-canonical ensemble treatment. Re-expressing the order parameters in terms of isospin 
density, we obtain (up to f/(l)/ phase), 

UA = 2F{m^nr)-^, {zuj'd) = -^(V^^)o ) (57) 



in agreement to leading order with previous result 

Thus, the appearance of the second order parameter ua is quite natural. Finally, we 
remark that the situation in the B phase is the mirror image of the above discussion. The 
new U{1)b breaking density is. 



iiu 



70^757-2^) = -AF^fieX (l - A') ^ cos(a) (58) 



C. 9 Dependence 

So far we have been mostly investigating the phase diagram in the {fis, f^i) plane at fixed 
6. In this section we would like to focus more on the 6 dependence, drawing the phase 
diagram in the {6, /i) plane. This trivial exercise leads to rather interesting consequences, 
namely, the 6 dependence at fixed /i becomes non-analytic. We further characterize the 
phase diagram in terms of the (GG) correlator and the topological susceptibility x- Finally, 
we confirm our calculations by checking the validity of low energy theorems. 

To simplify the discussion we shall take fij = and focus on 6 dependence in the Normal 
and Baryon phases. The situation in the Isospin phase is again just the mirror image, as 
can be explicitly checked. 
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We begin by considering 6 dependence at /i = 0. The story is exactly the same as in the 
well-studied case Nc = 3,Nf = 2. The vacuum energy density J^{9) is, 



^(^,/x = 0) = -AF'miie) (59) 



where, 



KiO) ^ (60) 

m{6) = ^ ((m„ + rridf cos'{d/2) + (m„ - rridf sm\d/2)) ^ (61) 
By differentiating J^{d) we can compute correlation functions of GG, 

dT ,.g^GG, 

^"^ - X=- /^^3;(T^^(x)^^(0))eonn (63) 



9^2-^ / ^ 327r2 ^ ' 327r' 



At /i = we find, 



/.g'^GG ImuTrid . ,7 ,\ 

(^^).^o = sm(^) 



The expressions flfi4j) reflect the well-known strong ^ dependence in the region m„ ~ = 
rriq, 6 ^ 71. Let's introduce the asymmetry parameter. 



\mu - md\ 
rriu + rud 



(65) 



and assume e <^ 1. 

The CP odd order parameter (iGG) (see Fig.Qa) starts out at when = and increases 
smoothly with 6, reaching its maximum just before 6' = tt at, 

(^^).^.---^(^^)o (66) 

Afterwards, the order parameter (iGG) experiences a steep crossover, dropping to its mini- 
mum of, 

^^ g^GG ^ 

The crossover occurs in the region — 7r| ~ e and hence, the topological susceptibility x has 
a sharp peak around 6* = tt of width A9 ~ e and height, x(^)/lx(0)l = (see Fig-IDb)- 

Such behaviour of the CP odd order parameter (iCG) strongly suggests that for rriu = 
rud, spontaneous breaking of CP symmetry occurs at 6' = vr. This situation, known as 
Dashen's phenomenon, has been extensively studied in A^^c = 3 QCD with Nf = 3 and 
Nf = 2j^ H, M, III- For Nf = 3 with ^ mu,md it is believed that spontaneous 



(^^)^-^- - +^(^^)o (67) 
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FIG. 1: e dependence in Nc = 2, Nf = 2 QCD at = 0, e = 0.01. a) The CP odd order parameter 
{i ^^Ff ). See eq. for precise normalization, b) Topological susceptibility x- 



CP breaking occurs at 6* = tt for |m„ — rridlms < rriumd. For Nf = 2, CP violation occurs 
at ^ = TT, rriu = and possibly in a small window of \mu — J^j-dl 7^ 0,11]. 

However, it is important to note that Dashen's phenomenon is not under complete theo- 
retical control in our effective Lagrangian lj^UI) . Indeed, for a moment, we fix m„ = m^. Then, 
for general 6', the mass term explicitly breaks the symmetry of the effective Lagrangian (j^Uj) 
from S'f/(4) to Sp{A). However, for 6 = tt, the mass term in the effective Lagrangian van- 
ishes, restoring the symmetry to SU{A) and giving rise to apparently massless goldstones: 
m^(6' = tt) = 0. Yet, no such symmetry restoration occurs in the fundamental micro- 
scopic QCD Lagrangian at 6* = vr. This contradiction is resolved by including higher order 
(quadratic) mass terms in the effective Lagrangian, which would explicitly break S'f/(4) even 
at ^ = TT [lO| . It is precisely these terms, which control the physics of Dashen's phenomenon, 
and which are not included in the present work. 

We do not wish to consider such higher order mass terms in this paper. For any fixed 
7^ these terms can be neglected by considering sufficiently small rUg. If the higher 
order mass terms are largely saturated by a third quark of mass mu,d <^ ^qcd, then 

we require, 

\mu - rridl i^u,d fnl{9 = 0) 

> — ^2 (68) 



V 



This condition is, indeed, realized in the true physical world. If, on the other hand, the 
higher order terms are controlled by a light rj' (as motivated by 00), we consider, 

\mu-md\ m{^i))^ mljO = 0) 

m^ + md ~ M^, ^^^^ 

Of course, by imposing restrictions ()68|) . (jU^ we automatically exclude the regions of param- 
eter space where Dashen's transition is realized, and we may discuss only the quantitatively 
steep crossover in the Normal phase. However, we shall see in a moment that by considering 
the system at finite /i, the rapid changes in the vicinity of ~ tt observed in the Normal 
phase will be washed out. 

Let us now turn on finite fiB- Once conditions (jUHj) . are met, all the results of 
previous sections hold for any 6. In particular, the transition to the Baryon phase occurs 
at /i = mT,{6) (see Fig. 12)). As explained above, we can consider arbitrarily small ratio 
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FIG. 2: Phase diagram of Nc = 2,N} = 2 QCD as a function of /i^ and 9. Here, e = ^^"^^ = 0.01 
A rapid crossover occurs in the Normal phase at ^ = vr, which is conjectured to become a first 
order phase transition, when m„ = m^. 



m^iO = n)/m'^{9 = 0) = e as long as — > 0. Thus, for /is < mj^^O = n) the Normal phase 
is realized for all 6, while for fiB > TnT^{9 = 0) we are entirely in the Baryon phase. Finally, 
if we fix /i^ with rri-^lO = tt) < /i^ < m-^lO = 0) and vary 9 from to 27r we encounter 
two phase transitions: from Normal to Baryon phase and then back to Normal. Thus, 
the 6 dependence becomes non-analytic in this region! Since the N to B phase transition 
is second order, we expect the topological susceptibility, x to be discontinuous across the 
phase boundary. The transitions between B and N phases occur at 6' = 6'c and 6 = 27i — 6c, 
with the critical 9c given by mT,{6c) = /i. 

In the Baryon phase, the free energy density reads, 

m = -2FVI (l + ^) (70) 

Clearly, the 6 dependence in the superfluid phase is different from that in the Normal phase 
()59|) . This is most clearly seen by computing, 

,.g^GG muTTid ,T 2 . 



X = -li^(^^)o'-^W (71) 



We have to remember that expressions (j71|) hold for all 6 only once we are entirely in the 
Baryon phase: /x^ > 171^,(9 = 0). On the other hand, if m^(6' = n) < hb < mT,{0 = 0), then 
we use expression ()64|) . for 6 where the Normal phase is realized, and expression (f7T|) . for 6 
where the Baryon phase exists. Focusing for a moment on /i^ > m,r(^ = 0), we see that the 
9 dependence is very smooth: there is no sign of rapid crossover in (iGG) near 9 = n and the 
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large peak in the susceptibility x disappears. Moreover, as hb increases, the 9 dependence is 
suppressed, as expected. This smooth 6 dependence at 6* ~ vr in the superfiuid phase should 
be contrasted with sharp behavior in the Normal phase discussed above, see Figs^a,b. 

Now we would like to understand, how the strong 6 dependence at /i = gets smoothed 
out as the chemical potential ^ib increases. For, < /i^ < m^(6' = n), 6 dependence is the 
same as at /i = 0. The key region is m^(^^ = n) < fis < rn.„{9 = 0), where at fixed fis, the 




FIG. 3: a) 6 dependence of the topological susceptibility X) i^i = 2, Nf = 2 QCD at fi = 
1.02 m7r(^ = vr) (broken curve) together with x at ^ = (unbroken curve). Here, e = 0.01 b) the 
same with = 1.3mT^{9 = vr); c) the same with /u = 0.7 mT^(9 = 0); d) the same with ^ = mT^{9 = 0) 



Baryon phase is realized for 9c < 9 < 2tt — 9c and the Normal phase is realized otherwise. 
Let us investigate the behaviour of the topological susceptibility x this region. As we 
expected, x is discontinuous across the phase transition, 

X{0t) -X{0c) _mumdml{0){ilj^)l . 2 



|X(0)| 6AF^fi% 



sm\9c) (72) 



As the chemical potential increases slightly past 171^^(9 = tt), a narrow region of Baryon 
phase appears around 9 = ir, deep inside the crossover region shown on Fig.^ In terms 
of susceptibility X; this affects only the very top of the peak of x(^) by introducing small 
discontinuities at 6* = 6'c and 9 = 2n — 9c (see Fig.Ela). As /i further increases, the range of 
9 where the B phase exists starts growing. This is accompanied by growth in discontinuities 
of X at the transition points. Eventually, for m^(6' = vr) < ^ m.,,{9 = 0), the original 
peak in x associated with the CP crossover, is entirely replaced by discontinuities associated 
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with the second order Normal to Baryon phase transition (see Fig.Olb). Once this occurs, 
the magnitude of the jump xij^t) ~ xi^^c) starts to decrease (Fig.Elc), until finally at /i = 
m.„{6 = 0), X becomes continuous again and we are entirely in the Baryon phase (Fig. Eld). 
The washout of the sharp 9 dependence near 6 = ir has been realized! 

The most exciting result of this section is that for ^ rrid the "Dashen's crossover" 
first splits into two second order Normal to superfiuid phase transitions and for fi > m,r(0) 
gets entirely washed out. We would like now to provide some speculations regarding the 
deg enerate case ruu = md, 9 = tt. This point might be of importance for lattice fermionsj^, 
15|, as it is equivalent to a theory where one quark mass is negative and 9 parameter is not 
explicitly present. In principle, it is possible to analyze this situation rigourously by going 
to higher order Chiral Perturbation Theory. However, since the algebra becomes rather 
involved even at /i = 0, we confine ourselves to a conjecture based on the above results and 
common 9 = n lore. For Nc =^i^Nf = 2, at /i = 0, we expect that Dashen's phenomenon 
will occur along the same lines[lO| as for N^. = 3. Spontaneous CP violation will happen at 
9 = Tc, however, no continuous symmetries will be broken and Goldstones will have a small, 
but finite mass 171^^(9 = vr) > 0. At finite fi, we expect a line of first order phase transitions 
at ^ = TT to extend to /i = m.„{9 = vr), where it sphts into two second order phase transition 
lines (see Fig.@]). We remark that in the Baryon phase, P-parity is still spontaneously broken 
at = TT, while in the Isospin phase, P-parity is broken at 6' = 0, but not a 9 = tt. 



<iGG> > 



<iGG> < 
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FIG. 4: Conjectured form of the Phase diagram of Nc = 2, Nj = 2 QCD for = m^. Solid line 
indicates a first order phase transition, while dashed Unes indicate second order phase transitions. 
The region near 6 = tt is not to scale. 



V. WARD IDENTITIES 



In this section we check the validity of Ward Identities (WI) i. El El Q at nonzero 
fi, 9. We anticipate that WI must remain untouched when external parameters such as /x, 
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9 or temperature T are introduced. Indeed, the anomaly (chiral and/or conformal) is a 
short distance {UV) phenomenon, which is not affected by medium effects (density /i 7^ 0, 
9 and/or temperature T). This fact was imphcitly used when we constructed the effective 
Lagrangian (j2Up . However, we are in a position to calculate each term entering the WI 
explicitly. Therefore, the check of the WI is a nontrivial test of self consistency of our 
results. 

The first identity that we consider, relates the two CP odd order parameters, 

(^1^) = ]^(^^75M^) (73) 

The identity (ffH|) reflects the well known fact that there is no 9 dependence when m — > 0. 
By consulting Table 1 and eqs. (jMj) . (fTTj) . we can explicitly check that our results satisfy 
the identity ()73p both in Normal and superfiuid phases. 
The next WI we would like to discuss is, 

/({inn a^CC 1 
"^'"^^^ W^"^^ |^(0))co„n = -^m^) + 0{M') (74) 

The 0{M'^) term in (j74p is usually dropped in the chiral limit, SU{Nf)y symmetric limit at 
9 = 0, assuming the resolution of the U{1) problem when flavor singlet 1]' is a heavy state. 
Indeed, in this case Table 1 and eqs. (jMj) .(f7T |) imply that the WI flTIj) holds both in Normal 
and superfiuid phases. An important remark is that both sides of (fTHj) and (ffij) depend on 
/i in a very nontrivial way. Nevertheless, the identities are preserved as expected. 

Now, we would like to see what happens with ()74|) when we relax the requirement of the 
SU{Nf)Y symmetric limit and also consider 9^0. In this case it is important to keep the 
0(M2) term, 

OiM^) = J d^xiTtfj'y.M^ix) ^75MV'(0))eo„„ (75) 
We begin in the Normal phase and evaluate, 

X--|WMV-) = -iM^(^*>. (76) 

The above result implies that singularities develop in the O(M^) term of eq. ()74|1 . due 
to r^'/goldstone mixing, which occurs for m„ 7^ m^. However, the singularities disappear 
when niu = rrid, so that the O(M^) term can be neglected in the chiral, SU{Nf)v limit, as 
long as we are sufficiently far from 9 = ir. This is the physically expected result. 

What happens with (fTH) in the superfiuid phase, when 9 ^ 07 We can immediately 
see that the O(M^) term can no longer be neglected even when m„ = m^. As has been 
discussed above, the Normal to superfiuid transition is second order, so that the topological 
susceptibility x generically experiences a jump across the phase boundary (f7^. while the 
chiral condensate (ipMip) (Table 1) is continuous. Thus, for the WI (f7i|) to hold, the O(M^) 
term must jump across the phase boundary, accounting for the discontinuity in x and, thus, 
contributing on the same footing as (ipMip) to the righthand side of (j74p . It is not surprising 
that the O(M^) term becomes important. Indeed, from Table 1 and (f7T|) we see that both 
X and {i/jMi/j) are of order M^/ /i^ in the superfiuid phase. The fact that x becomes O(M^) 
rather than 0{M) is part of smoothing out of 9 dependence in the superfiuid phase. The 
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contact between 0{M) dependence in the Normal phase and O(M^) dependence in the 
superfluid phase is provided by the fact that /i^ ~ 0(M) at the Normal to superfiuid phase 
transition. Thus, all correlators in eq. (fTIj) develop fi~^ singularities in the superfiuid phase, 
which are due to the modification of the goldstone spectrum (j2D). 
Finally, from dH}, to leading order in M, in the superfiuid phase, 

j dS{Ti,^^Mij{x) ^75M7/'(0)),o„„ = -J^^^K + K- 2m„mdCos(^))(V^V')o (77) 
which vanishes only if ttIu = met, 6' = 0. 



VI. GLUON CONDENSATE 

Having determined the 6 and dependence of different condensates and densities con- 
taining the quark degrees of freedom (Tables 1, 2), one can wonder if similar results can be 
derived for the gluon condensate (G^,^), which describes the gluon degrees of freedom. As 
is known, the gluon condensate represents the vacuum energy of the ground state in the 
limit ruq = 0, = and plays a crucial role in such models as the MIT Bag model, where a 
phenomenological "bag constant" B describes the non- perturbative vacuum energy of the 
system. The question we want to answer: how will the gluon condensate (G^^,) (bag constant 
B) depend on /i, 9 if the system is placed into dense matter? This question is relevant for a 
number of different studies such as the equation of state in the interior of neutron stars, see 
e.g.[l3, or stability of dense strangeletsjT^]- Of course, it is difficult to answer this question 
in full 3 color QCD at finite /i^, however, the answer can be easily obtained in 2 color QCD 
for /i <^ ^QCD, which is the subject of the present work. We limit ourselves to considering 
only the Normal and Baryon phases, the results in the Isospin phase, as always, are obtained 
by replacing /x^ ^ /i/. We work in Minkowski space in this section. 

We start from the equation for the conformal anomaly. 



e>; = -^^G%G^>^^ + ijMi: (78) 



where we have taken the standard 1 loop expression for the /3 function and b = yA'c — |A^/ = 
6 for Nc = Nf = 2. As ususal, a perturbative constant is subtracted in expression (f75|). 
For massless quarks and in the absence of chemical potential, eq. (I78|l implies that the QCD 
vacuum carries a negative non-perturbative vacuum energy due to the gluon condensate. 

Now, we can use the effective Lagrangian (j^n|) to calculate the change in the trace of the 
energy-momentum tensor (9^) due to a finite chemical potential /ib -C Aqcd- The energy 
density e and pressure p are obtained from the free energy density J-', 

e = + fisnB (79) 
p = -J^ (80) 

Therefore, the conformal anomaly implies, 



Gl,G^nt.,m,e-{^,Gl^G^^'"'), = -4(^(/i,m,^) -^o)-/iBriB(/i,m,e) + (V^Mz^), 



(81) 
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Here, the subscript on an expectation value means that it is evaluated at jj, = m = 0, 6 = 0. 
The good news is that we have already calculated all quantities on the righthand side of 
eq. (|HT|) - see expressions (jSUj) . (fTOj) and Tables 1, 2. Thus, in the Normal phase we obtain. 



i^G^C^n.^m^e - (^G;,G^-)o = -3mi9){,p^), (82) 



When ^ = 0, (I82p reduces to the standard result 11411 . which was derived in a different manner 



As expected, (G^^) does not depend on fi in the Normal phase. The Baryon phase is more 



exciting. 



{^,G'',.G>^n,,^,e - {^,G%G^no = (l + . (83) 

It is instructive to represent the same formula in a somewhat different way, 

, bg' GlGn,,m,e - {^,GlG^n,=o,n.,e = AF\^,l - mUO)) (l - 2^^) , (84) 



which makes contact with the fact that in the Normal phase, when /x^ < 771^^(6), the gluon 
condensate does not vary with fiB- However, for /i^ > ^n{(^), the dependence of the gluon 
condensate (G^j^) on fiB in the Baryon phase becomes rather interesting. The condensate 
decreases with for m.,, < jjLB < 2^/^m^ and increases afterwards. The qualitative differ- 
ence in the behaviour of the gluon condensate for jis ~ 'm-n and for <ti fiB ^qcd 
can be explained as follows. Right after the Normal to Baryon phase transition occurs, the 
baryon density is small and our system can be understood as a weakly interacting gas 
of diquarks. The pressure of such a gas is negligible compared to the energy density, which 
comes mostly from diquark rest mass. Thus, (Op increases with ub and, according to the 
anomaly equation (f7H|l . (G^^,) decreases. A similar decrease in (GL) with baryon density is 
expected to occur in "dilute" nuclear matter (see [3] and review Jlil|). On the other hand, 
for fiB 3> m^, energy density is approximately equal to pressure, and both are mostly due 
to self-interactions of the diquark condensate. Luckily, the effective Chiral Lagrangian ()20|) 
gives us control over these self-interactions as long as /i^ ^ ^qcd- Such control is largely 
absent in corresponding calculations of (G^^) in nuclear matter. As Ae ~ Ap, the trace 
(Op decreases and the gluon condensate increases with baryon density. Such behaviour of 
(G^j,) is quite unusual, as finite baryon density, on general grounds, is expected to suppress 
the gluons. 

At small baryon density, we can also use a slight variant of the above method for calcula- 
tion of the gluon condensate, originally developed in the context of nuclear matter. As long 
as our Bose-condensate of diquarks is dilute, we may neglect interactions between diquarks, 
and approximate the change in (G^^,) as just the expectation value of G^j, in each diquark 
times their number. 

Here q~ denotes a diquark state relativistically normalized to {p)\(l^ {p')) = 
2Ep{2'n')^S^{p — p'), giving rise to the factor g;^- in ()85|) . It remains to calculate the matrix 
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element, (g ). This can be done by sandwiching the anomaly equation (f7H|) between 

two diquark states. As {q~\Qa\(l~) = 2?ti^, 



2ml = -{q-\^G';,,Gnq-) + (g"|V^M^|g-) (86) 



We are used to the fact that the goldstone mass comes entirely from the symmetry breaking 
term t/jMi/j, so we might, naively, expect from eq. ()86p that vanishes in a diquark. 
However, the diquark mass, to first order in M is given by, 

ml = {q-\ijM^\q-) (87) 

therefore, 

{'l-\7;^G:^Gnq-) = -ml (88) 



so that and ipMip contribute equally to the goldstone mass in eq. (jHE}- Now from 
eq. (IHSD, 

{-^^G'^^uG'"'"') tM,m,e - {^^G"^uG^''"')o,mfi = --nsfn^, (89) 

This is in agreement with our full result ()83|) . (j84|) to leading order in n^. 

Finally, we note that by differentiating (j82|l . (jHBI) with respect to 6 we can obtain corre- 
lation functions of with GG, in Normal and superfiuid phases. 



VII. CONCLUSION. SPECULATIONS. 

We conclude this paper with the following speculative remarks: 

1. Naively, one could say that we studied in the present paper a purely academic 
question by considering = 2 rather than the realized in nature QCD with Nc = 3. 
We should comment on this as follows. First, for A^^ = 2, the so-called, diquarks become 
well-defined gauge invariant objects. However, diquarks, as has been argued in a number of 
papers, (see, e.g. recent papers on the subject fiO. I21I li^. I23I li^. lis^) may play an important 
role in 3-color QCD dynamics. If the passage from SU{2)coior to SU{'i)coior does not lead 
to dramatic disturbances of these diquarks, these predictions based on SU{2)coior remain 
qualitatively valid in real QCD! Arguments supporting the conjecture of smoothness of the 
transition SU{2)coior to SU{?>)coior are presented in 2J|. We also note that there is some 



similarity between the proposal of [25] and the present work to study the diquark dynamics. 
In the proposal 25| the idea is to introduce a color source to study the diquark dynamics, 
while in our paper with Nc = 2 the diquarks automatically become gauge invariant objects, 
and no source is required to study them. 

2. It has been suggested that the 6'-induced CP odd state can be formed in heavy ion 
collisions at RHIC, see original papers[26] and a recent review |27l|. Our analysis could be 
quite relevant for the study of the decays of a GP odd configuration, if it is formed. 

3. It has been known for quite some time that violation of parity and CP parity (which 
is the case when 6' 7^ 0) may completely change the phase structure of a theory. Some lattice 
calculations, for example, suggest that the behavior of the system could be very nontrivial 
when 6* 7^ 0, and some singular behavior and even phase transitions are expected, see e.g.(28l|. 
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In an environment where C, P and CP are strongly violated, the interaction of quarks and 
anti-quarks is not identical, as it is usually assumed, but rather, could be very different. 
Under such conditions the QCD phase transition in the early universe could have a much 
more complicated history than it is typically assumed. In particular, one can imagine that 
some very nontrivial objects, such as Witten's nuggets |29|. which behave as dark matter 
particles, can be formed. Moreover, due to the differences in interactions between quarks 
and anti-quarks in the presence of 6, local separation of baryon charges may take place, and 
chunks of quarks or anti-quarks in condensed color superconducting phase may form during 
the QCD phase transition, serving as dark matter [30|. This scenario is based on the idea 
that while the universe is globally symmetric, the anti-baryon charge can be stored in chunks 
of dense color superconducting antimatter. In this case, instead of baryogenesis, one should 
discuss the separation of baryon charges. Such a global picture of our universe, is definitely, 
not in contradiction with the present observational constraints |3fl| . Rather, it may give a 
natural explanation for some global parameters such as ^^da/ /Rb |3Q| , or even can naturally 
explain the famous 511KeV line from the bulge of our galaxy j3l|. 

Typical relaxation time for 9 is much larger than Aqcd^ therefore, one can neglect the 
dynamics of 6 for studying the possible phases for each given 6. This was the main reason 
for us to study 6 dependence of the QCD phase diagram. We find in the present analysis 
that, indeed, two phase transitions of the second order will take place when 6 relaxes from 
^ = 27r to ^ = 0. These phase transitions will occur under very generic conditions when /i/ 
is smaller than m^, but of the same order of magnitude as m^. The physical consequences 
of these phase transitions are still to be explored. 

4. Aside from these far reaching speculations, we would like to mention here a few much 
more terrestrial consequences of the present study, which may have some impact on the 
lattice simulations. First of all, 2 color QCD is a nice laboratory to study a variety of 
different very deep problems of gauge theories at nonzero temperature and density, see e.g. 
(3^ 33L 33 . 3^. New elements, which were not studied previously and which are the subject 



of the present work, are related to 9 dependence of different condensates. There are a few 
interesting observations which deserve to be mentioned here: 

a) At ^ = TT, when the determinant of the Dirac operator is real, the superfluid phase is 
realized at much lower critical chemical potential fic than at 6* = 0. In the limit m„ = m^, 

we expect, ^"1^"^! ~ (a^^) ^ ■ gives a unique chance to study the superfluid phase on the 
lattice at a much smaller /i than would normally be required. We hope that our conjecture 
for the disappearance of Dashen's transition in the superfluid phase can be explicitly tested 
on the lattice. 

b) Knowledge of 6 dependence of different condensates allows one to calculate the topological 
susceptibility and other interesting correlation functions as a function of /i. Corresponding 
Ward Identities at nonzero /i can be tested on the lattice. 

c) Physics of gluon degrees of freedom and fi dependence of the gluon condensate can also 
be tested on the lattice. The behavior of the gluon condensate as a function of fi is very 
nontrivial, as has been explained in the text. Nevertheless, our prediction is robust in a 
sense that it is based exclusively on the chiral dynamics and no additional assumptions have 
been made to derive the corresponding expression. 

Finally, we should emphasize that all results presented above are valid only for very small 
chemical potentials /i_B,/i/ -C Aqcd when the chiral effective theory is justified. For larger 
chemical potentials we expect a transition to a deconfined phase at f^sifJ^i) — '^■^qcd tSy]. 
We should also add that all results presented above can be easily generalized to Nc = 3 
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QCD with /ij ^ 0, /iB = ollj. 
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APPENDIX A: PARAMETRIZATION OF THE VACUUM MANIFOLD 

The purpose of this appendix is to clarify some global aspects of the manifold of goldstone 
fields in iVc = 2, Nf = 2 QCD. Once the vacuum manifold is parameterized, we show that 
to first order in quark mass, all effects of the 6 parameter can be incorporated into a common 
real quark mass. 

We begin with the assumption that the chiral symmetry breaking pattern of A^^c = 2, Nf = 
2 QCD is, exactly, SU{4:) 5p(4), so that the vacuum manifold X = SU{A)/ Sp{4:). We 
represent the manifold as, 

X = {UJ^cU^, U e SU{4)} (Al) 

with Ec given by eq. (jS]). Note that X G A, where ^ is the set of all 4 x 4 antisymmetric, 
unitary matrices of determinant 1. The original work||2] had implicitly assumed that X = A. 
As we shall show, this is almost, but not quite true. In fact, A = X (j iX, i.e ^ is a disjoint 
union of two pieces, both of which are homeomorphic to SU{4:)/ Sp{4:). 

Even though such technical details do not affect the analysis of they become im- 
portant, once the ^—parameter is introduced. In particular, if we minimized the effective 
Lagrangian ()2()|1 over S G we would obtain very different 6 dependence. In fact, the the- 
ory (j2Up with S G iX corresponds to the theory with S G A" with the redefinition 6 ^ d + n. 
As long as we represent our vacuum manifold by any one of the two pieces X or iX, we 
obtain the same physics, if we define 6 appropriately. However, if we enlarge the vacuum 
manifold to contain both pieces, the physics changes: we obtain cusps in 9 dependence at 
6 = instead of Dashen's phenomenon at = tt. Since we find no evidence of an addi- 
tional spontaneously broken discrete symmetry that would connect the two disjoint pieces 
of A, we insist on our original assumption that the vacuum manifold is SU{4)/ Sp{4) = X. 

Now, let us demonstrate the above claims. We begin by writing any E G ^ as. 

Here, a, b are complex numbers and we have used the fact that S is antisymmetric. The 
requirement, det S = 1, implies, 

{det{C) + ahf = 1 (A3) 

We call, A± = {T, ^ A, det{C) + ab = ±1}. Then A is the disjoint union, A = ^+ U^^. 
We shall show, A+ = X, A- = iX. We begin with the observation that X C A+. Indeed, 
SU{4:) is connected, so X is connected. But, Sc G X, and detiC) + ab = 1 for Sc. Therefore, 
det{C) + ah=l for any S G A", so A* C A+. 
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Now, take E G The condition = 1 implies, 

C^t + |a|2 = 1 (A4) 
CtC+ = 1 (A5) 
aVsC^aa = feC^ (A6) 

The remaining condition for S G ^+ is, 

(iet(C) +ab = l (A7) 

Equation ()A4|) imphes \a\ < 1, and C can be written (non- uniquely) in the form. 



C= y^l- \a\^e'^^S, S G SU{2) (A8) 

Substituting (jMl) into (jMI), produces, |a| = We write, a = \a\e''^% b = |6|e*'^\ If a = 0, 
the condition ()A6|1 is satisfied automatically, and ()A7|1 implies, e'^^'^" = 1. If |a| = 1, the 
condition (|X6| is again satisfied, and (jX7|) gives, e*('^''+'^'') = 1. Finally, if < a < 1, (|X6)) 
becomes, 

e'^-^^as^^as = e'^^'^+^''^ (A9) 
But St/(2) is pseudo-real: aa^^as = S^, so e^'"^- = e*("^"+<^''). Substitution into gives 

We note that if e'^'^" = —1, we can always reabsorb the negative sign into the definition of 
S G SU (2). Thus, S G ^+ if and only if it may be written as. 

We now show that any E of form ()A11|) is in X. We let, = y^l - |a|2 + i|a|. Define 
matrices, Ui, U2, U^,U G SU{A) as, 

f/i = exp fzi -rOl (A12) 



2 \ zo"2 

U2= ( 1 ) (A13) 

[/a = exp(i(/.„S/2) (A14) 

(A15) 

= f/sf/sf/i (A16) 

One can now check that UT^JJ'^ = S. This concludes the proof of the fact that X = A+. It 
is now trivial to show that A- = iA+ = iX. 

The most practically useful result of the above discussion is the explicit parametrization 
(jAllll of the vacuum manifold X. For instance, this parametrization allows one to prove that 
the local minimum of the effective Lagrangian ()20|) at finite fiB and fij, originally constructed 
in j^, is, actually, global. Moreover, we can now use the form ()A11|) to study the topology 
of the vacuum manifold X. Indeed, the matrix S G SU{2) can be uniquely written as 
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S = Xo + iXi(Ti, XqXq + X,Xi = 1. Also, e^'f'- = X4 + 2X5, X| + X | = 1. Clearly, A" is in 
one to one correspondence with the set of vectors in M^, a/1 — |ap(Xo, Xi, X2, X3, 0, 0) + 
|a|(0, 0, 0, 0, X4, X5), < |a| < 1. But this is just a parametrization of S^. Hence, X = 
SU{4)/Sp{4) ^ S\ 

Finally, let us discuss the SU{A) transformation (j^H). Since, UqBUI = B, UqIUI = I, 
the kinetic part of the effective Lagrangian (|2(Jj) remains invariant. Therefore, we have to 
discuss only the mass term: 

£^ = -gReTiiMi:) = -gReTiiMUotU^) (A17) 

where S G A" and, therefore, can be written in the form ()A11|) . as, 

\ -C-' a* (T2 / 

Expanding, 

J^m = ^9 (e~^^/2Tr(e*""H^A C}) + e'^/^TT{e-''"'^{M, C*})) (A19) 

At this point, we again use the pseudo-reality of SU{2), C* = a2Ca2, obtaining after some 
algebra. 

Cm = 2gm{e)TT{C) = -gm{e)TT{MoT,) (A20) 

Thus, to first order in M, all 6 and m„, dependence can be incorporated into the quark 
mass m{9) ((221) • 
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